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Abstract 

Are Markovian master equations for quantum Brownian motion independent of 
model assumptions used in the derivation and, thus, universal? With the aim of 
answering this question, we use a random band-matrix model for the system-bath 
interaction to derive Markovian master equations for the time evolution of one- 
dimensional quantum systems weakly coupled to a heat bath. We study in detail 
two simple systems, the harmonic oscillator and the two-level system. Our results 
are in complete agreement with those of earlier models, like the Caldeira-Legget 
model and, in the large-band limit, with the Agarwal equations (both with and 
without rotating-wave approximation). This proves the universality of these master 
equations. 
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1 Introduction 



The description of the interaction of an open quantum system with its envi- 
ronment is an important problem in quantum physics. If the environment is 
modeled as a heat bath, the interaction will lead to relaxation and dissipa- 
tion processes in the quantum system, and to an irreversible approach toward 
equilibrium jjj]. During the last decades, various models of this type have been 
introduced in different branches of physics and chemistry. Notable examples 
are the Redfield theory in nuclear magnetic resonance ||, the Oppenheim- 
Romero-Rochin model in condensed-phase chemical physics ||, the phase- 
space approach of Agarwal in quantum optics JH|, and the influence func- 
tional method used by Caldeira and Leggett in condensed-matter physics 0. 
The Markovian master equations obtained in these approaches have been re- 
cently compared in Ref. M. 
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The standard model for quantum dissipation (the Caldeira-Leggett model) 
consists of a quantum system coupled to an infinite set of harmonic oscillators. 
Caldeira and Leggett have shown that it is always possible to treat the bath 
as an ensemble of independent oscillators provided the system-bath coupling 
is weak. These authors also assumed that the coupling is linear in both, the 
position coordinate of the quantum system and the bath variables. For the 
quantum system, this choice follows from the requirement that in the classical 
limit, the friction force should be linear in the velocity. For the bath, the 
choice was primarily made for computational convenience. To the best of our 
knowledge, there is no compelling argument for choosing the interaction term 
linear in the bath coordinates. 

It is expected, of course, that the relaxation process described by a master 
equation should be insensitive to the detailed form of the system-bath inter- 
action. The work reported in the present paper aims at proving this statement. 
We do so with the help of an alternative model for the interaction. We use an 
ensemble of random matrices. The ensemble encompasses all forms of system- 
bath interaction which are linear in the position coordinate of the quantum 
system, and which respect fundamental symmetries of the problem like time- 
reversal invariance. The ensemble is characterized by a few parameters which 
establish the relevant time scales. The Markovian master equations derived 
in this fashion are then valid for all possible forms of the interaction between 
quantum system and heat bath, except for a set of measure zero. We find that 
in the the high-temperature limit, the Markovian master equations derived 
by Caldeira and Leggett and others are independent of both the specific struc- 
ture of the bath and of the specific form of the system-bath interaction. This 
establishes the universality referred to in the title of our paper. 

Random-matrix theory (RMT) was originally introduced by Wigner to de- 
scribe spectral fluctuations of quantum many-body systems such as nuclei 
and has since been applied successfully in a wide range of other fields such 
as quantum chaos and disordered mesoscopic systems [|J. To the best of our 
knowledge, a random-matrix approach to relaxation has been first developed 
in nuclear physics in the context of deeply inelastic heavy-ion collisions ||. 
Since then, related models have been used to describe relaxation of a non- 
degenerate two-level system [[TTJ, dissipation in complex quantum systems 
P"T1 , |T2|1 and, more recently, the dynamics of a simple quantum system in a 
complex environment ||13|| . 

A second motivation for our work relates to the use of random-matrix models 
in closed quantum systems with many degrees of freedom. In such systems, 
only few degrees of freedom usually command physical interest. We refer to 
such degrees of freedom (to the remainder) as to the collective (the remaining) 
degree(s) of freedom, respectively. A case in point is nuclear fission. Here, 
interest is focussed on the shape degree of freedom, and little attention is 
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usually paid to the intrinsic degrees of freedom of the fissioning nucleus. In 
cases like this, the dynamical behavior of the remaining degrees of freedom 
depends, however, strongly on that of the collective degree of freedom and, 
therefore, cannot be modeled as a heat bath. The success of RMT in self-bound 
many-body quantum systems || then suggests that we model the remaining 
degrees of freedom in terms of a suitable random-matrix model. We have in 
mind an ensemble of random matrices which depends parametrically on the 
collective degree (s) of freedom. Such an approach has been taken in the papers 
by Bulgac et al. Before investigating the consequences of such an idea, 
it is necessary to study the limiting case where the environment can indeed 
be modeled as a heat bath, and to ask whether in this case, the heat bath 
can be replaced by a suitable random-matrix model. Our proof of universality 
answers this question affirmatively. 

The paper is organized as follows. In Section we present the model for sys- 
tem plus bath with emphasis on the statistical properties of the interaction, 
a random band matrix. An approximate form for the second moment (shown 
later to be equivalent to the rotating-wave approximation) is introduced. In 
Section ||| we derive the Markovian master equation for the averaged density 
operator of the system. The range of validity of this equation is discussed in 
Section f|. In Sections |5] and [| we apply our results to the damped harmonic 
oscillator and the dissipative two-level system, respectively. We obtain a gener- 
alized band-width dependent fluctuation-dissipation relation. We show that 
in the large band-width limit we recover the Agarwal equations (with and 
without rotating-wave approximation). Conclusions are drawn in Section [7[ 



2 The model 



We study the properties of a small quantum system S coupled via a random 
interaction to a large environment, considered as a heat bath. The Hamiltonian 
of the composite system is given by 

H = H s ®^s + ^b®H b + Q®V = H + W (1) 

where Hg describes the system S (for example, a particle moving in a potential 
or a spin degree of freedom), Hb describes the bath B (the actual form of Hb 
is not specified) and W = Q <8> V the system-bath interaction. We denote by 
\n) (| a)) the eigenstates of the system (bath) Hamiltonian with eigenvalues 
E n (e Q , respectively), 

Hg\n) = E n \n) H B \a) = e a \a) . (2) 

The product states \na) form a complete set for the composite system. The 
operator Q acts on the system S, and V is a Gaussian random band matrix 
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acting on the bath. The first two moments of V are given by 



V a b = VabVcd = (5 ac 5 bd + 5 ad 5 bc )V ab 2 (3) 

where the matrix V ab respects time-reversal symmetry and has non-zero el- 
ements only in a band of width A along the diagonal. More specifically, we 
adopt a form first given in Ref. 0. This paper also contains a detailed jus- 
tification of the form of Eq. (01). This form has been widely used later, cf. 



Ref [12,13 
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Vab = A [p(e a )p(£ 6 )p e" 5a*- . (4) 

Here A is the strength of the coupling, p(e) the density of states of the bath, 
and A the band width. For W, this implies 



(I) WT = W2 n W p J = (5 ac 5 bd + 5 ad 5 bc )Q mn Q pq V ab 2 . (5) 
The form of the Hamiltonian ([!]) is a generalization of the Hamiltonian consid- 



ered in Ref. ||10|| . There it was motivated with the remark that in relaxation 
problems, the process is frequently found to be insensitive to the details of 
the interaction. One may therefore construct an ensemble of interactions and 
calculate the average of the observable over this ensemble. We show below 
that this is indeed the case. 

In Eq. (|5p, only the part of the interaction acting on the bath is modeled as a 
random matrix. This is physically sensible since only the bath is supposed to 
be a complex system. We observe that as a consequence, the variance (I) of W 
has the inconvenient feature of being not symmetric in the variables of both 
the system and the bath. Therefore, we also consider a symmetrized form (II) 
of W where the entire interaction behaves as a random matrix. This form may 
be thought of as an approximation to the full form (I). We show later that 
form (II) leads to the rotating-wave approximation. It is given by 



(II) WT = W$ n W% = (5 ac 6 bd 5 mp 5 nq + S ad 6 bc 6 mq 6 np )\Q mn \ 2 V ab 2 .(6) 

Let p be the density operator for system plus bath. The von Neuman equation 
for p reads 

p(t) = U(t)p(W\t) (7) 

where the time-evolution operator U (t) = e~ lHt obeys Dyson's equation 

t 

U(t) = U (t)-i f dhU Q (t - hjWUih) (8) 



and where U (t) = e lHot denotes the free time-evolution operator. (We put 
h = 1 throughout). We define the reduced operator for the system S by 
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Ps = t r B[p] where the trace is taken over the bath states. We assume that the 
interaction is turned on at t = and that at that initial time S and B are not 
correlated. Then 



p(0) = p s (0) ® Pb(0) (9) 

is the product of the initial density operators ps(0) and Pb(0) for system and 
bath, respectively. We suppose further that at all times t > the bath is in 
thermal equilibrium with temperature T, 

p B = i$>-^|a>(a| (10) 



where Z is the canonical partition function and (3 = (kT) 1 . Expression flHl| ) 
can be shown (see Appendix [A]) to be equivalent to 

Pb = \a*)(a*\ and p(e) = p^e 136 , (11) 

where the state \a*) is defined by the temperature T. 



3 Derivation of the master equation 



In this Section we derive a Markovian master equation for ps(t)- The equa- 
tion applies provided the coupling between system and bath is weak. More 
precisely, we use the following assumptions. 

i) The time t obeys the inequalities £a "C t <^ tp. Here t& = 1/A is the 
duration time of a single action of the interaction and tp is the Poincare 
recurrence time of the system. This condition is always needed to describe 
a relaxation process in terms of a transport equation. 

ii) For all states |n) and |6), the band width A has to satisfy the inequalities 
w,7< A< kT. Here 7 is the relaxation constant and u denotes the mean 
level spacing of the system S. For the harmonic oscillator, 7 is defined in 
Eq. (p8|). For other systems, an analogous definition applies. Condition ii) 
requires weak coupling between bath and system and ensures the validity of 
the Markov approximation. It also requires the temperature T to be larger 
than a minimum temperature kT m = A and may, therefore, also be seen as 
defining a semiclassical approximation. 

These assumptions are discussed in more detail in Section [|. 

Because of the stochastic nature of the interaction W, the time-evolution 
operator @ and, consequently, the density operator ([?]) are themselves random 
variables. We have to calculate their mean values. The averaging procedure 
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consists in expanding U(t) in Eq. (H) in powers of W (Born series), averaging 
term by term, and finally summing up the whole series, see Appendix [B]. We 
always work in the limit in which the dimension N of the bath matrices tends 
to infinity. We consistently omit terms of order iV _1 and smaller. 

We illustrate the procedure by calculating the transition probability per unit 
time. We suppose that S is initially in some eigenstate \m) with ps(0) = 
\m)(m\, and we ask for the probability to find the system in another state \n) 
at a later time t. We have 

Pn(t) = (n\p s (t)\n} = J2(nb\p(t)\nb} 

b 

= J2(nb\U{t)\ma)(ma\U\t)\nb) . (12) 

b 

Expanding U(t) and W(t) up to first order in W and taking the average 
(denoted hereafter by a bar), we obtain 

b \ &n + £b — &m — £a ) 

For t ^> £a, the factor 4(sm 2 |xt)/x 2 is sharply peaked at x = and may be 
approximated by 27rt5(x). This yields 

Pjt) — 2irt ^2 | (nb\ W\ma) | 2 8{E n + e b - E m - e a ) . (14) 
b 

The transition probability per unit time is defined as 



W nm = = 2nJ2\(nb\W\ma)\ 2 5(E n + e b - E m - e a ) . (15) 

t b 

This is Fermi's Golden Rule. 

We turn to the derivation of the master equation. The averaged density oper- 
ator p(t, t') for system plus bath obeys the equation 

p(t,t') = U(t)p(Op\t') 



+ Jdrj dr' U(t-r)Wp(T,T')WU\t' - r) . (16) 



The averaged time-evolution operator U (t) obeys 



U(t) = U (t) + fdhj dt 2 U (t - t^WUoih - t 2 ) WU{t 2 ) . (17) 
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To obtain an evolution equation for the averaged reduced density operator p s 
we have to (i) solve Eq. (|T7|), (ii) substitute the resulting U{t) into Eq. ( |16|) 



(iii) take the trace over the bath degrees of freedom and (iv) differentiate with 
respect to time. We shall see that it is not always possible to perform step 
(i). This is the case, in particular, for the unsymmetrized variance of Eq. (|5|). 
Nevertheless, a master equation can still be derived in the weak coupling limit. 
This is done in the next two Subsections. 



3.1 Case (II) 



We begin with the simpler case, i.e., with the approximate form (II). In this 
case, the averaged time-evolution operator U(t) is diagonal and reads (see 
Appendix g Eq. ( pTj )) 



U nb (t) = e -i{E n+ e b )t-^t 

The decay width Y nb is given by Eq. ( [B.12|) 



T nb = 27r^2\(nb\W\n 1 b 1 }\ 2 5(E n + e b - E ni - e bl ) . (19) 

nifei 

We substitute U(t) from Eq. ([17]) into Eq. (|16|) and arrive at 

^2(nb\p(t,f)\nb) = Y,e- t{En+£b){t - t ' ) -^ (t+t 'Hnb\p(0)\nb) 
b b 
t t' 



+ J dr J d r ' e (-^"+ £ >)-¥)('-) 



o o 



X 

bn-i b\ 



\{nb\W\n l b l )\ 2 {n 1 b 1 \p(r,r')\n 1 b 1 )e^ En+Eb) - * ^''^ . (20) 

bnibi 

We take the time derivatives and get 

(| + |y) X>%(M')K> = -J2T nb (nb\p(t,t')\nb) 

o 

+ ]T |(n6|W|ni6i)| 2 / dt 1 e^ E " +£ >' )+E ^ )tl (m&i|p(* + t u OK&i) 

bnibi 



+ KnblWln^}] 2 J ^eH^+^+xW^^I^^ + ^lnj^) . (21) 



bni&i 

We have put t 1 = r — t and t[ = r' — t' . 
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The right-hand side of Eq. ([H]) is easily interpreted. The first term (the "loss 
term") corresponds to transitions which deplete the state n whereas the two 
last terms ("gain terms") correspond to transitions which feed the state n. 
The gain terms are non-local in time and, therefore, involve memory effects. 
In the limit of weak coupling, however, the process becomes Markovian and 
the memory effects play no role |L4[] . To see this, we note that for short times 
(to zeroth order in W), Eq. flUf ) reduces to 

p(t,t') = U o (t)p(0)Uo i (t') . (22) 

We use this form in the gain terms and accordingly approximate (n 1 6 1 |p(t + 
ti,t')\nxbx) by e^^i+^i^ni&ilpfttOK&i) and (nih\p(t,1f + ti)|ni&i) by 
(n 1 b 1 \p(t,t')\nibi)e i( - En ^ +£b ^ t ' 1 . Setting t = t r , we obtain 



jJ2(nb\p(t)\nb) = -J2T nb (nb\p(t)\nb) 

b b 



Kn&IWlni&OlVW^i&i) 

bn\b\ 





X 



J dt 1 e ( * iiEn+£b ~ En i~ £b i )+E: % k)tl + J dtie i ~ i{En+£b ~ Eni ~ £b i )+L % k)tl 



-t 



.(23) 



Condition ii) implies r n & <C E n + e b — E ni — e bl . Thus, we neglect T n b in the 
integrands. This amounts to replacing U(t) by Uo(t). Moreover, for sufficiently 
large times we have 

t 

J dtie (i(E n +e b -E ni -e H ) tl ^ 2 n5{E n + E b - E n - S bl ) . (24) 
-i 

Hence, Eq. ( p3|) reduces to 



j t Y.{nb\p{t)\nb) 



2tt J2 \(nb\W\n 1 b 1 )\ 2 5(E n + e b - E n - s bl )(nb\p(t)\nb) 

briib\ 



-2tt]T Kn6|W|ni6 1 )r<J(S n +e 6 -S ni -e 6l )(7n6 1 |p(t)|7i 1 6 1 ) . (25) 

bnibi 

We have used Eq. (|19|) for T nb . Using Eq. ([TJ) we finally obtain a master 
equation of the Pauli type, 



dP (t) 

— If 1 = E X> a &i|?(f)l"i&i> - E W nin J2(nb\p(t)\nb) 
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= Y,W nni P ni (t) -Y, W mnPn(t) . (26) 
ni ni 

The transition probabilities are given by the Golden Rule expressions fll5D , 



W nni =2irJ2\(nb\W\n 1 b 1 )\ 2 5(E n + e b - E n - s bl ) 



W nin = 2ir'£\{nb\W\n 1 b 1 }\ 2 5{E n +e b -E n -e bl ) . (27) 

It may seem curious that in the evaluation of the gain term, it is necessary 
to invoke the Markov approximation while in the loss term, the limit of weak 
coupling apparently suffices. To explain this fact, we recall that we always work 
in the limit of infinite matrix dimension N. The loss term in the Pauli master 
equation is obtained from single-side Wick contractions, symbolically written 
as VV( ) : ( ) or ( ) : ( )VV, whereas the gain term is generated by Wick 
contractions V( ):( )V which connect matrix elements in both amplitudes. 
Selection among the first type of Wick contractions is affected by both, the 
limit N — > 00 and the weak-coupling limit. We exemplify this statement by th e 
terms of fo urth order, given by the three Wick contractions WW, WW, 
and WW. The lines indicate which pairs of Vs are Wick-contracted. Of 
the three, the last is neglected because N — > 00 and the second, because of 
weak coupling. In contradistinction, the form of the gain terms in Eq. ( |2T| ) is 
determined entirely by the limit iV — > 00. Hence, an additional step is needed 
to implement the weak-coupling limit. 



3.2 Case (I) 



For the non-symmetric form (j5|) of the variance of the interaction, the averaged 
time-evolution operator is not diagonal in energy representation. Instead, we 
have 

(nb\U{t)\n'b') = 6 w (nb\U(t)\n'b) . (28) 

Therefore, the matrix elements (nb\U(t)\n'b') cannot be given in closed form 
(see Appendix 0). However, the time derivative of these quantities can be 
obtained explicitly in the limit of weak coupling. This suffices to obtain the 
master equation. Aside from this difference, the derivation proceeds in com- 
plete analogy to that given in the previous Subsection. In analogy to Eq. fl2"0| ) 
we obtain 



J2{nb\p{t,t')\n'b) = ]T (n6|?7(f)|rii6)(n 1 6|p(0)|n 2 6)(n 2 6|?7 t (Ok'6) 

b bn\ri2 
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X 



t t' 

J dr J dr'J2Yl ( n b\U(t - r)|r2i6)Q nin2 Q n3n4 Vb fe2 2 

bn\ » l 2 fl 2 



o o n3 „ 4 



x(n 2 b2\p(r,T f )\n 3 b 2 )(n 4 b\u\t' - r')\n'b) . (29) 



We take the double time derivative, use Eqs. (p.25| ) and ( P-26|) for the time 



derivative of the time-evolution operator, the Markov approximation in the 
gain terms, and the identity 

dre lXT = IP- + ii5(x) . (30) 
o 

We neglect the level shift due to the principal-value integral. All this yields 
the master equation 



^{n\p s (t)\n') = -i(n\[H s ,p s (t)} W) 



dt 

-n E WnZntnii^lPsit)^') ~ 7j E W nX'n 2ni ( n \PS (*) K) 

~ E W^ in2nl ( ni \p s (t)\ri2) + \ E ^Lr>i|Ps(*)K> • (31) 

ni"2 niri2 

We have defined the generalized transition probabilities 



Wnnxmraa ~~ 2?1 ~ E QnmQmm ^661 2 ^(-^n 2 + ^ni _ 
WnmWi = 2?r E Qmn'Qnim Hbi 2 <K-En 2 + £&~ E ni - e bli 



^inman' = 2?r E QnmQnan'Hb! 2 5(E„ 2 + £ 6l — En'- £ 6 ) . (32) 



4 Discussion: Time Scales 



With an eye on the derivation given in the previous Section, we discuss the 
various time scales appearing in our model. These time scales play an essential 
role in defining the range of validity of the master equation |fL4"l,|l5t ■ 



According to the statistical ansatz in Eq. (|j), the interaction V connects eigen- 
states of Hb within an energy interval ~ A. Thus, the band width A can be 
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visualized as the average amount of energy exchanged during a single action 
of V, and t& = 1/ A can be interpreted as the duration time of a single action 
of V (i.e., the time needed to transfer the energy A). A statistical description 
in terms of a master equation can be valid only for times t^$>t&. 

Any dynamical process in a finite-sized system will return (close) to its initial 
state after a characteristic time, the Poincare recurrence time tp (in the exam- 
ple of a two-level system, the Poincare time corresponds to the Rabi period 
tp = 2tt/(E + — EJ), where (E + — EJ) is the energy difference between the 
two (perturbed) levels). When the bath is much larger than the system, the 
recurrence time is essentially determined by the mean level spacing D of the 
bath, tp ~ 1/D. Obviously, tp tends to infinity with the size of the bath. This 
is the condition of irreversibility. The inequality t -C tp must be fulfilled in 
order to have relaxation. 



The weak-coupling condition ii) requires that the relaxation constant 7 be 
much smaller than the amount A of energy transferred during a single action 
of the interaction, 7 <C A (see also Appendix [FJ). This condition has a simple 
interpretation in terms of the times that correspond to these energies. The 
relaxation time tp = I/7 must be much larger than the time £a = 1/A 
needed for a single action of V, t& <C t R . 



There are two time scales which determine the memory time of the heat bath 
[jT6,[T7l : The time t& (the inverse of the frequency cutoff of the bath in the 
Caldeira-Leggett model) and the time tp = 1/kT. The latter is purely quan- 
tum in origin. For high (low) temperature, thermal (quantum) fluctuations 
dominate and the memory time of the bath is given by t& (by tp, respec- 
tively). A crossover between thermal and quantum fluctuations occurs at the 
crossover temperature kT m = A. To garantee the validity of the Markov ap- 
proximation, the temperature must be much larger than the crossover tem- 
perature, A <C kT. This condition can be rephrased in terms of length scales: 
The range of the interaction must be much larger than the thermal de Broglie 
wavelength of the Brownian particle XdB = 1/ V^MkT '. 



5 First Application: Harmonic oscillator 



We illustrate our results for the case where the system 5* is a harmonic oscil- 
lator with mass M and frequency to. The Hamiltonian H$ is given by H$ = 
Tfiy + \Muj 2 x 2 ) and the energy spectrum reads E n = (n + 1/2) u, n = 1, 2, . . .. 
We assume a coupling linear in the position of the system, Q = x. We intro- 
duce the usual creation and annihilation operators a) and a. The elements of 
the matrix WVi n vanish unless \m — n\ = 1. 



11 



5.1 Case (II) 



The transition probabilities in Eq. (|27| ) are easily evaluated. We use Eq. (| 
For times t <C tp, we can replace the sum over b by an integral and obtain 



W nni =27vJ2\Qnn 1 \ 2 V bbl 2 5(E n + e b - E ni - e ni ) 

b 



2n\Q nni \ 2 A / de 



2 (gfc-^) 2 



:27r|Q nni |M e5^i-^e ^ 



(33) 



In the last line, we have used Eq. flTT|) for the density of states p(e). We use 
Q = x = J 2Mz( a + an d obtain for the only non-vanishing terms 



W nn _ x = W ne-*" W nn+1 = W (n + l)e* 

W n . ln = W ne*» W n+ln = W (n + l) e -f 1 

where 

... A 7T 



(34) 



(35) 



The master equation for the damped harmonic oscillator takes the form 



Pn{t) = Wnn-lPn-l(t) + W nn+1 P n+1 {t) - (W n - ln + W n+ln )P n {t) 

= W e^[(n + l)P n+1 (t)-nP n (t)} 
+ W e~^ [nP n -x(t) - (n + 

or, equivalent ly, 



(36) 



P n (t) = 2 1 (n th + 1) [(n + l)P n+ i(t) - 7iP„(t)] 

+ 2 7 n t?i [nP n ^{t) - (n + l)P n (t)] (37) 

where we have defined the relaxation constant 

2 7 = W (e2" - e~2 u ) , (38) 

and where n t h denotes the average number of quanta at temperature T, 

(39) 



1 



n t h 
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Eq. fl37|) coincides with the Lax-Louisell master equation for the harmonic 
oscillator in the energy representation and evaluated in the Rotating Wave 
Approximation (RWA) |L8]|§ 



= — iu>[a'a, p s ] + r ){2ap s a) — a)ap s — p s a^a) 

+ 2'jn t h(aps a ^ + a Vs a — aJdp s ~ Vs aa ^) ■ (40) 

The RWA amounts to neglecting rapidly oscillating terms in the interaction 
and is valid for weak damping, 7 <C uj. The coincidence between our result 
and the Lax-Louisell master equation does not extend to the time dependence 
of the non-diagonal elements (n\p s \n') with n 7^ n' of the reduced density 
operator. This is due to the form @ of the second moment which supposes 
that the entire interaction acts as a random matrix: In our approach, the gain 
term for these non-diagonal elements vanishes because of the Kronecker delta 
for the states of the system S appearing in condition (II). 

For the harmonic oscillator, Markovian master equations cannot simultane- 
ously fulfill the following three conditions: (i) The reduced density operator 
is positive definite for all times t > 0; (ii) for t — > 00, the reduced density 
operator attains thermodynamic equilibrium; (iii) in the classical limit, the 



equation is equivalent to a Langevin equation. This was shown in Ref. [19 



The master equation ( flOf) is of the Lindblad form [^OJ which guarantees the 
positivity of the reduced density operator. As a consequence, the quantum- 
classical correspondence with the Langevin equation is lost, however. 



5.2 Case (I) 



We first evaluate the generalized transition probabilities of Eq. (|32|). We pro- 
ceed in complete analogy to case (II). We find, for instance, 

^inUns = 2nA Q rmi Q nin J^- E ^e- ( -^ 2 ^ . (41) 
Using the explicit form of the matrix elements Q, we get 



Win-m-m = Wone^ W^ +ln+ln+2 = W ^(n + l)(n + 2) 

W^+m+m = W (n + l)e~^ W^_ ln _ ln _ 2 = W oy /n(n - l)e"^ (42) 

where Wq is given by Eq. (|35|). Proceeding analogously for the other transition 
probabilities and inserting the result into the master equation, we obtain 
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j t {n\p s {t)\n') = -iw{n - n'){n\p s {t)\n')) 
-jnth fy\jn(n - 1) (n - 2\p s \n') + (n + l)(n\p s \n') 
-j(n th + 1) (n(n\p s \ri) + ^{n + l)(n + 2) (n + 2|p s |n ')) 
-771*7, (Jn'(n' + 1) (n|p 5 |n' - 2) + (n' + l)(n\p s \n')J 
-j(n th + 1) ^'(nlp^n') + sjjn 1 + l)(n' + 2) (n|p 5 |n' + 2) 
+l n th (V nn'(n — l\p~ s \n — 1) + \Jn{n' + — l|p 5 |n' + 
+7(nth + 1) (y(n + l)n'(n + l|p s |n' - 1)) 
+7n i?i (^/rm/(n - l\p s \n' - 1) + \](n + l)n'(n + l\p s \n' - 1)^ 
+7(% + 1) \^Jn{n' + l)(n - l|p 5 |n'^ 

+2 7 (n t/ , + 1) ^/(n + l)(ra' + l)(n + l|p s |n' + 1)) (43) 

It is easy to check that this result coincides with the master equation in energy 
representation derived (without RWA) by Agarwal for a harmonic oscillator 
linearly coupled to a bath, 



dps 



iuj [aJa, p s ] 



dt 

— 7(0^ ap 5 — 2a~p s a) + ~p s a)a + a 2 ~p s — ap s a — a^p 5 a^ + ~p~s a ^ ) 
- inth (2[a\ [a,p s ]] + [a\ [a\p s ]] + [a, [a,p s ]}) . (44) 

The bath was modeled as an infinite set of harmonic oscillators, and the pro- 
jection operator technique was used ||. 

5.3 Fluctuation-Dissipation Theorem and High-Temperature Limit 
With the help of Eq. Q55D, Eq. (BST) for 7 takes the form 



A ir . ,P 

7 = e smn—u . (45) 

' Muj 2 v ; 

This equation expresses the fluctuation-dissipation relation. We see that a 
small band width A tends to exponentially decrease the damping coefficient. 
In view of our ansatz for the interaction between system and bath, this fact is 
not surprising. In the limit of large band width, A ^> u, Eq. reduces to 
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the Agarwal fluctuation-dissipation relation |4[] between the friction coefficient 
7 and the diffusion coefficient D, 

D Ssl , P 

7 = — — - e 2 sinh—uj . (46) 
M.u) 2, 

Here the diffusion constant D is determined by the strength of the coupling, 
D = A 7ie-^ . (47) 

In the high-temperature limit flu <C 1, Eq. (f|(]) reduces to the Einstein rela- 
tion 

D , ^ 

7 = . (48) 

1 2MkT K ' 



We also note that using a = (y^-x + i\J ^i^p) i n Eq. (fi"4D we find the 
Caldeira-Leggett master equation 



d Ps 
dt 



-iu [H s , p s ] - vy [x, {p, p s }] - D [x, [x, p s ] . (49) 



This equation does not have the Lindblad form. This implies that the positivity 
of the reduced density operator can be violated for certain initial states. In 



Ref. pll ) it was shown that positivity is guaranteed provided the dispersion 
a xx =< x 2 > — < x > 2 of the initial wave packet obeys the condition a xx > 



^dB 2 ■ 



6 Second Application: Two-level system 



We derive the master equation for a two-level system with upper (lower) level 
|+) (|— ), respectively). We introduce the Pauli spin matrices a x , a y and a z . 
The Hamiltonian of the system takes the form 

H s = ^^o^ (50) 

where u is the energy separation between the two levels. We write E n = 
^ooon, n = ±1 and take Q = a x . For the generalized transition probabilies we 
find 



r(l) 



W y _l + _=W e* UJ0 (51) 

and similar expressions for , and Here Wq is given by Eq. (|35|) 
with ijj replaced by u and M set to unity. Thus, 
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J t (+\Ps\+) = -^W e^(+\p s \+) - ~W e?"°(+\Ps\+) 

+\e-*" {-\p s \-) + \e--^{-\p s \-) 
= -2j(nth + l)(+|p s |+) + 2 1 n th (-\p s \-) (52) 

and, in a similar way, 

_2 7 n th (-|/5 s |-> + 2 7 (n tfc + l){+\Ps\+) , 
iwo(-\p s \+) ~ l(2n th + 1) [(-\p s \+) ~ (+\Ps\-)) , 
-iu (+\p s \-) ~ l^ n th + 1) [(+\Ps\-) ~ (-\Ps\+)\ ■ ( 53 ) 

Here n th is given by Eq. (|39D with u replaced by u Q . Simple manipulations 
with the spin matrices show that the master equation for the two-level system 
may be written as 



d l I- I 

d 
dt 
d / 



-) 



HPsl+> 



-) 



— p s = -2j(n t h + 1) a + o-_p s o- + cr_ + 2^n th cr+p 5 cr_ 

-iu o + oJp s o-0 + - j(2n t h + 1) [a + aJp s a^a + - a + p s a + ] 
+iuj a_a + p s (T + (j_ - "f(2n th + 1) [a_a + p s a + a_ - cr_p s a_] 
-2jn th cr_a + p s a-a + + 2^{n th + 1) a-p s a + (54) 

or, more simply, 



fcPs = —2 U o\. (7 *>Ps] ( 55 ) 

+ 7 n th (2(T + p 5 a_ + (T+p s (T + + CT-P 5 CT- - p S (T-(T+ - V-(T+p S ) 

+ ^{n th + l)(a + p s a + + oJp s o_ + 2a_p s cx + - cx + a_p s - p s a + a_) . 

We have defined the operators <r_ = |— )(+| and a + = |+)( — |. For n t h = 0, 
Eq. fl5"5| ) reduces to the Agarwal equation for spontaneous emission of a two- 
level atom ||. It may seem surprising that Eq. (|55l) applies for T = where the 
inequality A <C T is clearly violated and non-Markovian effects are present. In 
Ref. |22j it was shown that the high-temperature master equation is obeyed 
even at low temperatures provided the time t is larger than the memory time 
of the bath, t tg. This implies that condition ii) of Section (|3|) may be 
replaced by the less restrictive condition tg <t <C tp- 
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7 Conclusion 



Starting from a random band-matrix model for the system-bath interaction, 
we have derived a Markovian master equation for the average reduced density 
operator of a one-dimensional quantum system. We have assumed that the 
system-bath interaction is linear in the position coordinate of the quantum 
system, and we have considered two cases: (I) That part of the interaction 
which depends on the bath variables is a member of an ensemble of random 
matrices of proper symmetry; (II) the entire interaction is a member of an 
ensemble of random matrices of proper symmetry. The form of the master 
equation differs in both cases. The equation is valid in a domain of parameter 
values specified by the inequalities i) and ii) of Section |[ 

We have applied the master equation to two cases, the damped harmonic oscil- 
lator and the dissipative two-level system. For the damped harmonic oscillator 
and case (I), we have obtained the same equation as Argawal who considered 
a harmonic oscillator coupled to a heat bath which was modeled as an infinite 
set of harmonic oscillators. In the limit of high temperature, this equation 
coincides with the master equation of the Caldeira-Leggett model. For case 
(II), our master equation for the diagonal elements of the reduced density 
operator is identical to the Lax-Louisell master equation evaluated in the Ro- 
tating Wave Approximation (RWA). This is because in case (II) we impose 
conditions upon the interaction matrix elements of the position coordinate of 
the quantum system. These conditions are tantamount to the RWA. For the 
non-diagonal elements of the reduced density operator, these same conditions 
imply the vanishing of the gain terms. In this point our result differs from 
that obtained by Lax and Louisell. For the two-level system and case (I), our 
master equation reduces, at T = 0, to the Argawal equation for spontaneous 
emission of a two-level atom. 

We conclude that Markovian master equations for quantum Brownian motion 
derived in the weak coupling limit possess universal validity: These equations 
are independent of the specific microscopic model used for their derivation. 
This is not true of approximations like the Rotating Wave Approximation. 
Typically, such approximations violate certain invariance requirements (trans- 
lational invariance in the case of the RWA) and, thereby, lose universal validity. 



A Appendix 



The equivalence of Eqs. ([[(]) and ([□]) is closely related to the equivalence of 
the microcanonical and the canonical ensemble in the thermodynamic limit 
23fl . We thus consider a thermodynamical system in contact with a 



N 
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heat bath. The (canonical) partition function is given by 

OO CO 

Z{(3) = J dep{e)e- pe = J dee~ p£+lnp ^ . (A.l) 
o o 

Since e and In p(e) grow with N, the integral can be evaluated by a saddle- 
point approximation. Expanding the integrand up to second order around its 
maximum e*, we get 

Z{p)=^)e-^Jdee^^(e-er{^^) J- (A.2) 
The distribution in energy is a Gaussian centered at e* with a width 

Ae 



' d 2 lnp(e)' 



kT 2 C v . (A.3) 



We have used the fact that k In p[e) is the microcanonical entropy. For N — ► 
oo, Ae/e* ~ l/y/N becomes negligibly small and the Gaussian approaches a 
5-function. Hence, 

Pb(0) = ^E^KH - IO(a*l ( A -4) 

where |a*) is the eigenvector corresponding to the eigenvalue e*. Moreover, 
according to Eq. ( |A.2|) the density of states can be approximated locally by 



p(e*) = p e^ . (A.5) 

B Appendix 

B. 1 The average propagator 

The propagator K(t) is defined by 

K(t) =U{t)9{t) =e- mt 9(t) (B.l) 
and obeys Dyson's equation 

CO 

K{t) = K (t) - i J dtxK (t - t x )WK{tx) (B.2) 

— CO 

where K (t) = e- iHot 6{t) is the free propagator and 6(t) the unit step function. 
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To calculate the average of K(t), we use the energy representation and intro- 
duce the following pair of Fourier transforms 



OO OO 

G(E) = - [ dte lEt K(t) , Kit) = / dte- iEt G(E) . (B.3) 

— OO — OO 

The transforms of (|B.1| ) and ( |B.2| ) are then given by 
and 

OO 

G = G + G WG = G (WG y , (B.5) 

s=0 

respectively. 

To average G, we use Eqs. (|5],|6]) and the rule that the average over a product 
of Gaussian distributed W's equals the sum over all ways of Wick-contracting 
pairs of W's. We show below that in the limit of weak coupling, only con- 
tractions between adjacent pairs of W's have to be taken into account. We 
find 



G = ]T G (WG ) S = G + G WG WG + G WG WG WG WG + ■■■ 

= G 1 (B.6) 

1 - WG WG y ' 



B.l.l Case (II) 
Since 

(nb\G Q (E)\n'b') = 5 nn i5 w — — t-=— ——— = 5 nn ,8 w (G ) nb (B.7) 

& — {ti n + E b ) + IE 

and 



(nb\WG W\n'b') = Y^ W b ^{G Q ) n , H W h y' ( G ) n >v 



■ 5 nn /8bb'(Go) n b W^ 1 (G ) nibl , (B-8) 



the matrix elements of G become 
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{nb\G\n b ) 



E - (E n + e b ) - Y^mbx (G ) nibl 
= 6 ™' Sw (B.9) 

E — (E n + E b ) — Rnb(E) 

which shows that G is diagonal in the unperturbed energy basis. Here, 

R nb (E) = A nb (E) - i^jp- (B.10) 
where 



W nni2 

a,,,(/:) I'Y. 



T nb {E) = 2vr ^ W^~ 2 5{E - (E m + e bl )) . (B.ll) 



mbi 



For sufficiently high temperature, T nb (E) depends slowly on energy and can 
be approximated by 



r nb {E)~r nb {E n + e b )=T 



nb 



= 2nJ2W bb n 1 1 5(E n + e b -E ni -e bl ) . (B.12) 

For the averaged propagator, this yields 

G nb (E) = — ^- (B.13) 

1 ' E-(E n + e b )+i T -f y } 

where we have neglected the level shift A nb . Transforming back to time rep- 
resentation, we obtain 

K nb (t) = e -W+^)-¥*0(i) . (B.14) 

In the limit of weak coupling, contractions between non-adjacent pairs of 
W's are negligible. This can be shown by inspecting the contributions from 
the various products of Wick contractions. We consi der h ere the simplest 
and compa re p4| the imaginary parts of (nb\ n'b') and (nb\ 



case 



]Vg Wg Wg^v I n'b') (the real parts represent level shifts and do not con- 
tribute to the decay width). We are thus led to compare T nb (E) with 

Sni&i W b n b ™ l2 (Go)n b r nibl (E). These quantities are calculated as in Section || 
Assuming that |Q nni | 2 vanishes unless the states n x and n are close in energy, 
we find that in the high-temperature limit uo <C A <C kT, T nb (E) may be 
replaced by T nb (E n + e b ), and that in the weak-coupling limit 7 € A, the 
second term may be omitted in comparison with the first. 
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B.l. 2 Case(I) 

In this case, the propagator is not diagonal in the energy representation and 
we have instead 

(nb\G\n'b') = 5 w {nb\G\n'b) . (B.15) 
As a consequence, the series (|B.6|) is not geometric any more and cannot be 



summed easily. We thus have to find another method to evaluate the aver- 
age propagator. We actually calculate only the time derivative of the average 
propagator which is sufficient for the derivation of the master equation. 

Iterating the Dyson equation (P-5| ) once 

G = G + G W + G WG WG (B.16) 

and taking the average 



G = G + G W + GqWGqWG = G + G WG W G (B.17) 



yields the same result as in Eq. ( |B.6| ). We use this fact to average the time- 
evolution operator. In the interaction picture, we have 

.dUjft = w(t)U(t) U(0) = U(0) = 1 (B.18) 
at 

where 

U(t) = e iHot U(t) and W(t) = e iH otWe~ iHot . (B.19) 

This equation can also be written in integral form 

t 

U{t) = l-if dhWihpih) . (B.20) 





Inserting Eq. ( |B.2(J| ) into Eq. ( |B.18[ ) and averaging, we find 



d(nb\U(t)\nb) 
dt 

t 



dh (nblW&lmbJimhlWitjlnih) (n 2 b 2 \U{tx)\nb) . (B.21) 



We have W£?W%£ = S bb2 Q nni Q niri2 V bbl z , hence 
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d{nb\U(t)\n'b) _ r i(E n -E n2 )t t e i(E n2 +e b -E ni - £bl )r 
dt J ^ 

X QnniQniri2 Vbbi 2 (n 2 b\U(t - r)\n'b) (B.22) 

where we have set r = t — 1±. In the weak-coupling limit dU /dt is small and 

U(t — r) may be replaced by U (t). This approximation amounts to neglecting 
memory effects and is thus a Markov approximation. We obtain, 



d{nb ^ n ' b) =-W W^ inin2 ^-^ (n 2 b\U{t)\n'b) . (B.23) 
We have defined 



WW mn2 = 2nY J Qnn 1 Qn l n 2 V bbl 2 8{E n2 + e b - E n - e bl ) . (B.24) 
In the Schrodinger picture, 



i(E n + e b )(nb\U(t)\n'b) 

\ E W<$ inin2 (n 2 b\U(t)\n'b) . (B.25) 

nin.2 



i(E n , + e b )(nb\W(t)\n'b) 

\ E Wi% inin , (nb\W(t)\n 2 b) (B.26) 

niri2 

with 



W^ n , n2ni = 2nY / Qn 1 n'Qn2n 1 V bb MEn2+e b -E ni -e bl ) . (B.27) 
61 



B.2 The evolution equation for the average density operator 

The general density operator at times t, t' is given by 

p{t,t')9(t)9(t') = K{t)p{0)K\t') (B.28) 



d{nb\U(t)\n'b) 
di 



We have similarly 



d(nb\W(t)\n'b) 
di 
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where we have used the definition (|A.2 ) for the propagator. We take the double 
Fourier transform and obtain 

p{E, E') = G{E)p{ti)G\E') (B.29) 

where G(E) is given by Eq. ( P-5| ). Hence 



p(e, £i = EE G (WG y P (o)(GlwyGl . (b.so) 



=0 r=0 



In the weak-coupling limit, G appears sandwiched between two Ws which are 
contracted across p(0), and the general term in the expansion ( |B.30| ) attains 
the form 



GW- ■ -WGp(0)G^W- ■ -WG ] . (B.31) 
Thus, 



p = Gp{0)G f + G WGp{0)G j WG f 

+ GWG WGp ( o'jG^WG^WG^ + ■■■ (B.32) 

or 



p = Gp(0)G t + GWGJ)G ] WG ] . (B.33) 

Eq. ( |B.33| ), together with Eq. ( |B.6| ) for G, completely determines the dynamics 
of the system. Transforming back to the time representation, we can rewrite 
these equations as 

p(t,t')e(t)9(t') = K(t)p(0)K\t') 

oo oo 



J drj dr'K(t - r)Wp(r, T')e(r)e(r')WK\t' - r') (B.34) 



-oo — oo 



and 



oo oo 



K{t) = K (t) + J dh J dt 2 K {t - t^WKoih - t 2 )WK{t 2 ) . (B.35) 



-oo — oo 
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